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A hadronic cascade model based on resonances and strings is used to study mass dependence of 
relativistic nuclear collisions from p+Be to Au+Au at AGS energies (~ lOAGeV) systematically. 
Hadron transverse momentum and rapidity distributions obtained with both cascade calculations 
and Glauber type calculations are compared with experimental data to perform detailed discussion 
about the importance of rescattering among hadrons. We find good agreement with the experimental 
. . . data without any change of model parameters with the cascade model. It is found that rescattering 

is of importance both for the explanation of high transverse momentum tail and for the multiplicity 
. of produced particles. 
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I. INTRODUCTION 



There is now a big interest in studying strongly interacting matter at high density and/or temperature which is 
created in high energy nuclear collisions. Indeed, at high densities and/or temperatures, QCD predicts the chiral 
symmetry restoration and quark dcconfinement. How can we create such matter ? At present, high energy heavy ion 
collision is considered to be a unique way to create such dense and hot matter at laboratory. In order to find such a 
^— V , new form of nuclear matter, several heavy ion experiments have been and are being performed with Si(14.6A GeV/c) 

■ or Au(11.6A GeV/c) beam at BNL-AGS and with the O(200A GeV/c), S(200A GeV/c) or Pb(158A GeV/c) beam 
O ■ at CERN-SPS. 

On ■ Since high energy heavy ion collisions lead to a huge number of final states, many event generators have been 
proposed to explore these high energy nuclear collisions, with the aid of Monte-Carlo realization of complex processes. 
(~| In these event generators, there are mainly three categories of models. The models in the first category assume 
Glauber geometry for the treatment of AA collisions. For example, FRITIOF [|, LUCIAE |], VENUS §, HIJING [|, 
DPM 11, HI JET lH and LEXUS belong to this category. Final interaction among hadrons are included in VENUS, 
HIJET and LUCIAE. In these models, main quantum features during the multiple scattering are preserved within the 
^ ] eikonal approximation, and efficiently fast calculations are possible. However, these approaches are mainly designed 
for the extremely high energy collisions {y/s > lOAGeV). 

The models in the second category (parton cascade models), such as VNI § and ZPC |^, have been recently 
developed to implement the interaction among partons to study space-time evolution of partons produced in high 

■ energy nuclear collisions. These models have been originally designed to describe ultra- relativistic heavy ion collisions 
at collider energies, such as BNL-RHIC and CERN-LHC, and they have met some successes in describing heavy ion 
cohisions at CERN-SPS energies 

The third category of models is a transport model which is often referred to as 'hadronic cascade'. For example, 
RQMD |l|,|l|, QGSM ARC |l|l, ART (it), UrQMD Q and HSD (l| can be categorized here. They have been 
successfully used to describe many aspects of high energy heavy ion collisions in a wide range of incident energies. 
For the description of AA collisions in hadronic cascade models, the trajectories of all hadrons as well as resonances 
including produced particles are followed explicitly as a function of time. Nuclear collisions are modeled by the sum 
of independent hadron-hadron (hh) collisions without interferences. Two particles are made to collide if their closest 
distance is smaller than y^a{s)/TT, here a{s) represents the total cross section at the cm. energy ^/s. As a result of 
the hh collision, secondary particles will be produced according to the specific model with some formation time. One 
of the most distinct difference among these models may be in the way of implementing hadronic degrees of freedom. 
In RQMD and UrQMD, many established hadronic resonances are explicitly propagated in space-time, while ARC, 
ART and HSD do not include higher hadronic resonances. Although both modelings seem to give similar results if 
we see the final hadronic spectra inclusively, we expect thermodynamic quantities like pressure or temperature before 
freeze-out predicted by those models would be different from each other Another difference is the treatment of 
multiparticle production. String model is adapted in RQMD, QGSM, UrQMD and HSD, while in ARC and ART, 
final states are sampled according to the direct parameterization of the experimental data. The hadronic cascade 
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model based on the string phenomenology implies that some partonic degrees of freedom play some roles in reaction 
dynamics implicitly. In fact, the estimation of partonic degrees of freedom has been done recently within UrQMD [^ . 



ARC 16 1 has shown that 'pure' hadronic model can describe the data at AGS energies. At collider energies, however, 
explicit treatments of partonic degrees of freedom will be necessary. 

The main purpose of this work is to perform systematic analyses of collisions from pA to massive AA systems at 
AGS energies, for which high-quality systematic experimental data are available within the hadronic cascade 

model, JAMl.O, which has been developed recently based on resonances, strings and pQCD. 

The main features included in JAM are as follows. (1) At low energies, inelastic hh collisions are modeled by 
the resonance productions based on the idea from RQMD and UrQMD. (2) Above the resonance region, soft string 
excitation is implemented along the lines of the HIJING model (3) Multiple minijet production is also included in 
the same way as the HIJING model in which jet cross section and the number of jet are calculated using an eikonal 
formalism for pcrturbative QCD (pQCD) and hard parton-parton scatterings with initial and final state radiation 
are simulated using PYTHIA program. (4) Rescattering of hadrons which have original constituent quarks can 
occur with other hadrons assuming the additive quark cross section within a formation time. Since these features 
of the present hadronic cascade model, JAMl.O, enables us to explore heavy ion collisions in a wide energy range, 
from lOOA MeV to RHIC energies, in a unified way, it is a big challenge for us to make systematic analyses in these 
energies. In this paper, we focus on the mass dependence of the collision system at AGS energies. Other applications 
at higher energies are found elsewhere p^ . 

The outline of this paper is as follows. We will present a detailed description of cross sections and modeling of 
inelastic processes for hh collisions in section ^ because elementary hh processes are essential inputs for the hadronic 



cascade model. In Sec. Ill , we first study the transverse momentum distributions of protons, pions and kaons in p+Be, 
p+Al, p+Cu, p+Au, Si+A, Si+Cu and Si+Au collisions at the laboratory incident momentum of 14. 6A GeV/c. We 
discuss the role of rescattering by comparing the cascade model results with the Glauber type calculations. We then 
discuss the collision dynamics for truly heavy ion colliding system Au+Au collisions. The summary and outlook are 



given in Sec. IV 



II. MODEL DESCRIPTION 

In this section, we present the assumptions and parameters of our model together with the inclusive and the 
exclusive hh data including incident energy dependence. 

A. MAIN COMPONENTS OF THE MODEL 

The main components of our model are as follows. (1) Nuclear collision is assumed to be described by the sum 
of independent binary hh collisions. Each hh collision is realized by the closest distance approach. In this work, no 
mean field is included, therefore the trajectory of each hadron is straight in between two-body collisions, decays or 
absorptions. (2) The initial position of each nucleon is sampled by the parameterized distribution of nuclear density. 
Fermi motion of nucleons are assigned according to the local Fermi momentum. (3) All established hadronic states, 
including resonances, are explicitly included with explicit isospin states as well as their anti-particles. All of them 
can propagate in space-time. (4) The inelastic hh collisions produce resonances at low energies while at high energies 
( XL 4GeV in BB collisions 3GeV in MB collisions and >, 2GeV in MM collisions) color strings are formed and 
they decay into hadrons according to the Lund string model . Formation time is assigned to hadrons from string 
fragmentation. Formation point and time are determined by assuming yo-yo formation point. This choice gives the 
formation time of roughly 1 fm/c with string tension k = IGeV/fm. (5) Hadrons which have original constituent 
quarks can scatter with other hadrons assuming the additive quark cross section within a formation time. The 
importance of this quark(diquark)-hadron interaction for the description of baryon stopping at CERN/SPS energies 



was reported by Frankfurt group |12 19[]. (6) Pauli-blocking for the final nucleons in two-body collisions are also 
considered. (7) We do not include any medium effect such as string fusion to rope medium modified cross 

sections and in-medium mass shift. All results which will be presented in this paper are those obtained from the free 
cross sections and free masses as inputs. 
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B. Baryon-baryon interactions 

Let us start with the explanation of the resonance model for baryon-baryon (BB) collisions implemented in our 
model in detail. We assume that inelastic BB collisions are described by the resonance formations and their decays 
below CM. energy ^/s < 4GeV and at higher colliding energies, string formation and their fragmentation into hadrons 
are included based on a similar picture to that in the RQMD and the UrQMD model . The total and elastic 
pp and pn cross sections are well known. Fitted cross sections and experimental data are shown in Fig. 0. Inelastic 
cross sections are assumed to be filled up with the resonance formations up to -y/s =3-4GeV. At higher energies, 
the difference between experimental inelastic cross section and resonance formation cross sections are assigned to the 
string formation. 
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FIG. 1. The fitted total and elastic pp and pn cross sections which are used in the code together with measured data taken 
from particle data group |Q. CERN-HERA and COMPAS group parameterization are used at high energy 

The following non-strange baryonic resonance excitation channels are implemented for the nucleon-nucleon scatter- 
ing in our model: 

(1) NN -> 7VA(1232), (2) NN -> NN*, (3) NN -> A(1232)A(1232), 

{4)NN-^NA*, (5) NN N*A{1232), (6) ^ A(1232)A*, (7) ^ iV*iV*, 

{8)NN^N*A*, {9)NN^A*A*. 

Here N* and A* represent higher baryonic states below 2 GeV/c^. The pp and pn cross sections are calculated from 
each isospin components cr(/) (in some cases we ignore the interferences between different amplitudes): 

|2|/-.^j, u rM2 

/ 



a{hih2 ^ hahi) =}_^\C{hih2,I)\'\C{h3h4,I)r<j{I) (1) 



where C{hihj,I) is isospin Clebsch- Gordon coefficients. For N* and A* productions, the sum of production cross 
sections of several resonance species (A^(1440) ~ A^(1990) for N* and A(1600) ^ A(1950) for A*) are parameterized, 
and resonance species are chosen afterward (see below). The strength of each branches a{I) are determined from 
the exclusive pion production data |25| . Isospin 1 = 1 component for NN collisions can be extracted from the pp 
reactions. We assume that isospin 7 = components are determined from the pn reactions, then explicit form of cross 
sections in different isospin channels can be written down as follows, 

a{pp^pA+) = ^a{I =1), aipn-^nA+) = ^a{I = 1), (2) 

aipp^nA++) = ^a{I=l), a{pn ^ pA^) = ^<j{I = 1), (3) 

a{pp^pp*)=a{I ^l), a{pn^np*) = ^<j{I=l) + ^a{I = 0), (4) 
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a{pp ^ A+A+) = \a{I = 1), a{pn ^ A"A+) = ^<j{I = 1) + = 0), 
5 20 4 

a{pp ^ AOA++) - ^a{I ^ 1), a{pn ^ A-A++) - = 1) + \a{I = 0) . 

The functional form for the non-strange baryonic resonance formation cross sections is assumed to be 



d2 



(5) 
(6) 

(7) 



All parameters except one-A production cross section are listed in tables || and |l| for each isospin channel where all 
cross sections are given in mb and \/s^i^ denotes a threshold. One-A production cross section (t{NN iVA(1232)) 
is parameterized with the following functional form 



ai{NN A^A(1232)) 



0.0052840v/\/s/2. 0139999 - 1 28.0401(Vi/2.124 - 1) 



0.480085 



- 2.11477)2 + 0.01714052 ((Vs/2.06672) - 1)2 + 0.5764222 ' 
in order to ensure correct threshold behavior. Pionic fusion cross section {pp — > dn^) has been fitted as 

0.14648(Vi/2.024 - l)0-20807 0.12892(Vi/2.054 - l)"-08448 



a{pp — > dTr~^) 



((yi/2.13072) - 1)2 + 0.0424752 ((yi/2. 18138) - 1)2 + 0.0592072 ' 



(8) 



(9) 



In actual simulations, we have effectively included the cross section of the NN ird into the A production cross 
section for simplicity. Similarly the s-wave pion production channels NN NNtTs are simulated as the 7V(1440) 
production. 

TABLE 1. Resonances cross section parameters for 1=1, tts denotes s-wave pion production. 



Channel 




a 


b 


c 


d 




cri(AiV -> 


NN") 


24.94700 


2.48150 


2.63330 


0.425358 


2.162 


ai (NN 


A(1232)A(1232)) 


7.63181 


1.41140 


2.67784 


0.311722 


2.252 


ai (NN ~* 


NA") 


8.01615 


2.74161 


3.34503 


0.259703 


2.340 


ai {NN 


Ar*A(1232)) 


13.14580 


2.06775 


2.75682 


0.247810 


2.300 


ai {NN 


A(1232)A*) 


19.63220 


2.01946 


2.80619 


0.297073 


2.528 


ai iNN ^ 


N*N*) 


11.67320 


2.31682 


2.96359 


0.259223 


2.438 


0-1 (AA ^ 


AT* A*) 


2.99086 


2.29380 


3.54392 


0.090438 


2.666 


ai {NN ^ 


A* A*) 


35.13780 


2.25498 


3.14299 


0.215611 


2.804 


ai (NN 


NNtts) 


15.644100 


1.675220 


2.07706 


0.658047 


2.014 



TABLE II. Resonances cross section parameters for 1=0, tvs denotes s-wave pion production. 



Channel a bed V^th 

ao{NN NN*) 166.60600 2.10128 2.34635 0.284955 2l62 

aoiNN -> A(1232)A(1232)) 39.99770 1.83576 2.40348 0.288931 2.252 

aoiNN A(1232)A*) 56.32490 2.00679 2.71312 0.362132 2.528 

ao{NN N*N*) 2.14575 0.21662 3.40108 0.252889 2.438 

o-o(A^A ^ A* A*) 4.14197 1.67026 3.75133 0.476595 2.804 

aoiNN^NNns) 78.868103 0.746742 1.25223 0.404072 2.014 
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FIG. 2. The resonance production cross sections for pp (left panels) and pn (right panels) reactions as functions of cm. 
energies. In the upper panels, the total one-resonance (NR), double-resonance (RR), total resonance (NR-I-RR), total inelastic 
cross section and string formation cross sections are shown. In the middle and the lower panels, each one-resonance production 
branches and double-resonance production branches are plotted. 

In Fig. |2[ We show the contributions of non-strange baryonic resonance cross sections for different partial channels 
as functions of cm. energies. The upper panels of Fig. || show the one-resonance production cross section NN NR 
(solid lines) , two- resonance production cross section NN RR (dotted lines) , the sum of NR and RR cross section 
(long dashed lines) for pp (left panels) and pn (right panels) reactions. Total inelastic cross sections are filled up by 
the resonance productions up to about Ecm =4GeV, while at CERN/SPS energies, string excitation is dominated. 
The dot-dashed lines in the upper panels of Fig. express the string excitation cross sections for pp and pn. At AGS 
energies corresponding to the invariant mass y/s^ 5GeV, the contributions of the resonance productions and string 
productions are approximately the same in the first nucleon-nucleon collision in our parameterization. The collision 
spectrum in BB collisions, however, are spread in broad energy range for Au-|-Au collision as shown in Ref. Jigf , due 
to the high baryon density achieved at AGS energies. Low energy cross sections, therefore, is also important in order 
to treat the dynamics correctly at AGS energies from first chance NN collisions to the final hadronic gas stage. 

The cross section for the resonance productions may be written by 



dai2^34 (253 -fl) (254 + 1) 



PsilMfAimDAimDdimDdiml) , (10) 



dQ 647r2spi2 

where Si,i = 3,4 express the spin of the particles in the final state. Mass distribution function A{mi) for nuclcons is 
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just a (5-function, while that for resonances is given by the relativistic Breit-Wigner function 



2\ 1 '7iflr(m) 



^ J\f {m^ -ml)^ +ml.r{my ^^^^ 

where Af denotes the normahzation constant. In this paper, we use simply take Af = tt which is a value in the case 
of a constant width. The full width T{m) is a sum of all partial decay width TujMB) for resonance R into mesons 
M and baryons B which depends on the momentum of the decaying particle [ p^pS] ] : 

r^{MB) = TliMB)"!^ ( P^mf^Y^' (12) 



m \Pcms{mR)J 1 , n O / Pa^.irn) ^ ^^+^ 

where £ and Pcms (m) are the relative angular momentum and the relative momentum in the exit channel in their rest 
frame. 

The Monte Carlo procedure is as follows. First, final resonance types A(1232), N*oi A*are chosen using parame- 
terized cross sections and then we determine each resonance production channel according to the equation Eg. (p^ . 
To do this, we need to know the matrix element for all resonances. In the present model, we make a simple 

assumption that each resonance production cross sections can be selected according to the probability: 

P(i?„ R,) ^ {2S, + l){2Sj + ^)J J P^jA,{mj)Aj{'m])d{m:i)d{m]) . (13) 

Namely, the partial cross section for each resonance state is only governed by the final spins and mass integrals 
ignoring the resonance state dependence of the matrix element. Once types of resonances are chosen, we generate the 
resonance masses according to the distribution neglecting the mass dependence of the matrix elements in Eq. ( p^ ) 

P(to3, m4)dm3(irn4 ^ 4m3TO4P34A(m3)A(m4)(im3dm4 , (14) 

in the reaction 1 + 2 — > 3 + 4, where mass distribution function A{mi) should be replaced by the 5- function in the 
case of stable particles in the final state. 

In the case of the collisions involving resonance states in the ingoing channel, we use the approximation that the 
inelastic cross sections for resonance productions as well as the elastic cross section are the same as the nucleon- 
nuclcon cross sections at the same momentum in the cm. frame, except for de-excitation processes, NR —f NN and 



RR NN, whose cross sections are estimated by using the detailed balance described in the next subsection, II C 
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FIG. 3. The energy dependence of the exclusive pion production cross sections for proton-proton and proton-neutron 
interactions as a function of cm. energy. Sohd lines are the results obtained from the code. Data from Ref. |25| . 

Fig. H shows the energy dependence of the exclusive pion production cross sections in pp (up to five pion production) 
and pn (up to two pion production) reactions. We compare the results obtained from our simulation with the data [p5| . 
Overall agreement is achieved in these exclusive pion productions within a factor of two with the above simplification 
of the common matrix element in Eq. (|l3|). Smooth transition from the resonance picture to the string picture at 
Ecm = 3 ~ 4 is achieved since no irregularity of t he energy dependence is present in the calculated results. String 
excitation law will be described later in section [IE. In order to get more satisfactory fit, for example, we can improve 



the model by introducing different values for the matrix elements for different resonance channels. For example, in 
Ref. prf , the matrix elements are fitted to reproduce the pion production cross sections up to two-pion productions 
as well as rj production cross section assuming that they are independent of masses but dependent on species. 



C. Detailed balance 

In the processes of resonance absorption, we use a generalized detailed balance formula [p8|-p0t which takes the 
finite width of the resonance mass into account. The time- reversal invariance of the matrix element leads to the 
principle of detailed balance. If scattering particles are all stable, the formula is given by 
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rf«T34^12 i^Si + 1)(2S'2 + 1) PI2 dai2- 



»34 



dn (2S'3 + l)(2S'4 + l)p^4 dn 



(15) 



where Si denotes the spin of particle i and corresponds to the cm. momentum of the particles i and j. 
The differential cross section for the reaction (1,2) (3,4) for the stable particles may be written by 

rfcri2^34 |-Ml2^34p 1 ^ f f 2 2\,f 2\xi 2 2\it 2\ f^a\ 

M.^s {2S, + l){2S, + l) ^J J P^^Sip,-m,)dirn,)Sip,-mMm,) (16) 

where |A^i2^34p represents the spin-averaged matrix element. If the particles have a finite width, we should replace 
above 6 functions to the certain normalized mass distribution functions A{in). Using the |A^i2^34| — |A^34^i2|, we 
obtain 

rfg34^i2 ^ (2^1 + 1)(2^2 + 1) pfa dai2^34 1 . . 

dn (2^3 + l)(2^4 + l)P34 dn J J puA{ml)A{ml)d{ml)d{ml) ' ^ ' 

where we use the relativistic Breit-Wigner function Eq.(pT|) for mass distribution hmction A{m'^). The extra factor 
compared to the usual detailed balance formula increases the absorption cross section. It has been proved that this 
formula plays an essential role in order to understand the tt A^A dynamics [ps]-^ . 

In the case of one-A(1232) absorption cross section, we can write down the following formula: 



1 1 Pn 



,2n 7 2 



(^NA^NN' ^ 7;t~^ o-ww'^jvA / pA.(',o-}'2)A{'m)dm . (18) 

2 1 + SnN' pa \Jim,,+m^y J 

Pn and pA are the final nucleon-nucleon cm. momentum and the initial cm. momentum, respectively. The factor 
1/(1 + Snn') in Eq.(|T^) arises from the identical nature of the final states, and 1/2 comes from spins. 

There are some versions of the extended detailed balance formula, which are slightly different from each other. For 
example, Danielewicz and Bertch p8[ use the formula 

rfo'34^12 _ (25i + 1)(25'2 + l)pi2 ms rrn dai2^34 1 

dn ~ (25-3 + 1)(2S'4 + 1) P34 dn J J p34A{m'^^)A{m2)dm'^^dm'^^ ' ^ ' 

here denotes the pole mass of the resonance i, while Wolf, Gassing and Mosel use 

rffT34^12 ^ (2^1 + 1)(2^2 + 1) Ph dai2^34 1 



dn (25-3 + 1) (254 + l)p^4 dn J J A{ml)A{ml)dmldm: 



2 • 



(20) 



We have checked that these formulae give similar results to ours. Fig. |4| shows the comparisons between the different 
formulae of the cross sections for the reaction A++n pp. 
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FIG. 4. The cross section for A'^^n — > pp calculated using the different descriptions for the detailed balance as a function 
of ^/s — Ma, where Ma is the mass of ingoing A. The short-dashed lines correspond to the results of the formula which does 
not take the A width into account. The results of the cross sections obtained using the formula in Ref. ^] and are shown 
by dotted and full lines respectively. The long-dashed lines correspond to the results using Eq.(p^. 



D. Meson-Baryon, Meson-Meson Collisions 



We now turn to the explanation of meson-baryon {MB) and meson-meson [MM] collisions. We also use reso- 
nance/string excitation model for MB and MM collisions. 
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FIG. 5. Parameterization of the total and elastic ix^p FIG. 6. Parameterization of the total and elastic 7r~p 

cross sections. The data has been taken from |Q. Total and cross sections. The data has been taken from ]24[| . 
elastic cross sections are assumed to be the sum of s-channel 
and t-channel resonance and/or string formation processes. 
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Total cross section for ttN ingoing channel is assumed to be decomposed to 



(Ttot{sY^ = crsw(s) + (Tel{s) + Crg_g(s) + cr^_s(s), (21) 

where aei{s)^ (Jbw{s)^ crg_g(s), and cr^;_g(s) denote the f-channel elastic cross section, the s-channel resonance forma- 
tion cross section with the Breit-Wigner form, the s-channel and t-channel string formation cross sections, respectively. 
We neglect the t-channel resonance formation cross section at a energy range of y/s^ 2GeV. The t-channel elastic cross 
section (Jei{s) was determined so that the sum of the elastic component of the s-channel Breit-Wigner cross section 
aBw{s) and i-channel elastic cross section (Jei{s) reproduces the experimental elastic data for TriV interaction. Above 
the A(1232) region, t-channel elastic cross section becomes non-zero in our parameterization (Figs. ^ and P). String 
formation cross sections ((Jg_g(s) and (T^_g(s)) are determined to fill up the difference between experimental total 
cross section and (JBw{s)-\-(Tf,i{s). We calculate the resonance formation cross section (Tbw{s) using the Breit-Wigner 
formula |3^,|l^ (neglecting the interference between resonances) , 

^ ' Plm ^ i2SM + l)i2SB + l)iV^-mR)^+TRitot)y4 ^ ' 

The momentum dependent decay width Eq. (^2|) are also used for the calculation of decay width in Eq. (^2|). Sr, Sb 
and 5*^/ denote the spin of the resonance, the decaying baryon and meson respectively. The sum runs over resonances, 
R = A^(1440) ~ A^(1990) and A(1232) A(1950). Actual values for these parameters are taken from the Particle 
Data Group j2^ and adjusted within an experimental error bar to get reasonable fit for MB cross sections. The 
results of the fit are shown in Figs. ||,^ It has been shown that inclusion of resonances play an important role 
to study strangeness productions in AGS and SPS energies Q and {K~,K^) reactions In fact, strangeness 

production cross sections for resonance-iV ingoing channels are found to be much larger than that for ttN channel. 
This would be effective to explain the strangeness enhancement observed in heavy ion collisions within a rescattering 
scenario. 
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FIG. 7. Parameterization of the total and elastic K~p 
cross sections. The data has been taken from |2^]. The sum 
of the s-channel resonance productions and t-channel elastic, 
t-channel charge exchange and t-channel hyperon production 
cross sections can describe the data properly at low energies. 
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FIG. 8. Parameterization of the total and elastic K~n 
cross sections. Parameterizations of low energy part of the 
cross section where data is absent are addressed by Ref . . 
The data has been taken from [241. 



Since the KN interaction has some exoergic channels such as KN — > ttY, we need to include additional terms 



crtot{s) 



(Tbw{s) + <Tei{s) + CTch(s) -l- cr^y(s) -f crs-S(*) + '^t-S('^)' 



(23) 



where ach(s) and (Jt^y{s) denote t-channel charge exchange reaction and t-channel hyperon production cross sections 
which are also fixed by the requirement that the sum of i-channel contributions and Breit-Wigner contributions 
reproduce experimental data. Breit-Wigner formula enables us to calculate experimentally unmeasured cross sections 
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such as pN — > KK . For the calculation of crBw{s), we include hyperon resonances, R = A(1405) ~ A(2110) and 
I](1385) '-^ E(2030). The total and elastic cross sections for KN interactions used in JAM are plotted in Figs. and 
^ in comparison with experimental data ||2^ . 



K-p^7i:°A K-p^Tf^E- K-p^7i:°E° 




FIG. 9. The energy dependence of the exclusive hyperon production cross sections K' p tt' A, K^p 7r~S+, 
K~p TT K~n ix~YP , K~n — > tx^TP , and charge exchange cross section K~p K^n used in the model are 
shown by the full lines together with s-channel (dotted) contributions. Circles are data from Ref. |26| . 

The symbol (7.^y{s) in Eq.(p3|) is the sum of f-channel pion hyperon production cross sections KN — > ttY , y = A, S. 
In Fig. y, we plot the cross sections of hyperon productions and charge exchange cross sections as well as Breit-Wigner 
contributions fitted in our model. The cross section for the inverse processes such as t:Y — s- KN are calculated using 
the detailed balance formula. 
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FIG. 10. Parameterization of the total and elastic K'^p and K^n cross sections as well as charge exchange cross section in 
K'^n interaction. The data has been taken from [ p^ . 

KN ingoing channel cannot form any s-channel resonance due to their quark contents. Therefore the total cross 
section can be written within our model as follows, 

where cr|-_j^(s) is i-channel resonance formation cross section. Total, elastic and charge exchange cross sections 



used in our model are shown in Fig. |10| In the present version of JAM, only KN KA, KN K{892)N and 
KN i4r(892)A are explicitly fitted to experimental data |^ and fitted results are shown in Fig. |l^. 
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FIG. 11. Parameterization of the A, K* productions cross sections in K^p and K'^n interactions. The data from Ref. 
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In meson- meson scattering, we also apply the same picture as that in meson-baryon collisions: 



crtotis) = (TBw{s) + o-t-R(''^) + '^ei(s) + c^s-S(*) + '^t-S(*)- (25) 

The difference between experimental inelastic cross section and resonance cross sections at energies above resonance 
region for the meson-baryon and meson-meson collisions are attributed to the string formation cross section where 
l/^/s energy dependence of crg_s(s) is used |l3|. 

For the cross sections for which no experimental data are available, we calculate the total and clastic cross sections 
by using the additive quark model |^,|l2|,^ . 

2/3/ ATTV / 2/3 N /r,7N 

where a^N, express nucleon-nucleon total and elastic cross sections and Ui^Usi are the number of quarks and 

s-quarks contained in the hadron respectively. This expression works well above resonance region where the cross 
section becomes flat. 

For the t-channel resonance production cross sections f^.j^ls), we do not fit experimental data explicitly in this 
work except for NN reaction and one and two pion productions in KN reaction, because of the vast body of the 
possibilities for the final states. Instead, we simply determine the outgoing resonance types according to the spins 6*3, 
S'4 in the final state and phase space for the production of resonances -R3 and R^^ 

PiRs^Ri) cc {2S3 + l){2Si + l)p3iisf . (28) 

where ^34(5) denotes the cm. momentum in the final state. If the ingoing channel involves resonances, their ground 
state particles are also considered in the final state. Once the outgoing resonance types are determined, we generate 
masses according to the Breit-Wigner distribution. 

For the angular dependence in the processes of i-channel resonance production o'^_j^(s)j we use 

'-^^-eMbt), (29) 

and the slope parameter b for the energy range of ^/s > 2.17GeV is parameterized by 

6 = 2.5 + 0.71og(s/2), (30) 

with invariant mass squared s given in units of GeV^ . We use the same parameterization presented in Ref. |]3^ for the 
energy below y/s < 2.17GeV for the ^-channel resonance productions. The elastic angular distribution is also taken 
from Ref. || for < lOGeV and from PYTHIA || for > lOGeV. 



E. String formation and fragmentation 

At an energy range above y/s > 4 — 5GeV, the (isolated) resonance picture breaks down because width of the 
resonance becomes wider and the discrete levels get closer. The hadronic interactions at the energy range 4-5< ^/s <10- 
lOOGeV where it is characterized by the small transverse momentum transfer is called "soft process", and string 
phenomenological models are known to describe the data for such soft interaction well. The hadron-hadron collision 
leads to a string like excitation longitudinally. In actual description of the soft processes, we follow the prescription 
adopted in the HIJING model [Q, as described below. 

In the center of mass frame of two colliding hadrons, we introduce light-cone momenta defined by 

p+=E+p,, p-=E-p,. (31) 

Assuming that beam hadron 1 moves in the positive z-direction and target hadron 2 moves negative z-direction, the 
initial momenta of the both hadrons are 

Pi = {pt^Pi^^r), P2^ {pt,P2^^T) ■ (32) 
After exchanging the momentum {q^ tPt)j the momenta will change to 
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p[ = iil-x+)P+,x-p-,PT), p'^ = ix+P+,{l^x-)p-,-pT), (33) 

where P~^ = + = P^ = Pi + P2 = \^ (ii^ c.m. frame). The string masses will be 

^x~{l-x+)s-p^, ^x+{l~x-)s~p'^, (34) 

respectively. Minimum momentum fractions are x^^^^^ = /P^ and x^-^^ = p^ / P^ . For light-cone momentum 
transfer for the non-diffractive events, we use the same distribution as that in DPM g and HIJING |]: 



for baryons and 



^^""^^ (a;±2+cVs)i/4((i_a;±)2 + cVs)i/4 ' ^^^^ 



for mesons, where c — O.lGeV is a cutoff. For single-diffractive events, in order to reproduce experimentally observed 
mass distribution dlVP/M^, we use the distribution 

^(^*) = (a;±2+cVs)l/2- 

The strings are assumed to hadronize via quark-antiquark or diquark-antidiquark creation using Lund fragmentation 
model PYTHIA6.1 [ p2[ . Hadron formation points from a string fragmentation arc assumed to be given by the yo-yo 
formation point |^5|] which is defined by the first meeting point of created quarks. Yo-yo formation time is about 
Ifm/c assuming the string tension k = 1 GeV/fm at AGS energies. 

The same functional form as the HIJING model for the soft prp transfer at low pr < po is used 

/(Pt) = {(p| + c?)(p|+P^)(l + e(^'--^'°)/=^)}'' , (38) 

where ci = O.lGeV/c, po = 1.4GeV/c and C2 = 0.4GeV/c, to reproduce the high momentum tail of the particles at 
energies Eiab = 10 - 20GeV. 
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FIG. 12. The rapidity distributions for protons (circles), tt^ (triangles) and tt (squares) in pp collisions at 12GeV/c (upper 
panel) and 24GeV/c (lower panel). Histograms are the results obtained from our calculation. The data are from Ref. [M. 



14 



In Fig. y_2|, the calculated rapidity distributions of protons, positive and negative pions for proton-proton collisions 
at 12GeV/c and 24GeV/c are shown with the data from Ref. The proton stopping behavior and the pion yields 
are well described by the present model. Within our model, fast protons come from resonance decays and mid-rapidity 
protons from string fragmentation of Lund model (PYTHIA6.1 j2^) with the default parameter which determines 
the probability of diquark breaking at these ener gies . Anisotropic angular distribution in a resonance decay is taken 
into account assuming Gaussian pT distribution |13|] with a mean value of 0.35GeV/c^. 



p( 1 2GeV/c) + p p(24GeV/c) + p 




PT(GeV/c) pT(GeV/c) 



FIG. 13. The transverse momentum distributions for protons (squares), 7r+ (triangles) and vr (circles) in pp collisions at 
12GeV/c(left panel) and 24GeV/c(right panel). Histograms are the results obtained from our calculation. The data are from 
Ref. Ig. 

As reported in HIJING Q , an extra low px transfer to the constituent quarks is important to account for the high 
Pt tails of the pion and proton distributions at energies piab ~ 20GeV/c^. As shown in Fig. ^ the present model 
also reproduces the proton and pion transverse momentum distributions reasonably well at piab=12GeV/c^ and 24 
GeV/c2. 




E„,(GeV) E,„(GeV) (GeV) E„„ (GeV) 



FIG. 14. The energy dependence of the inclusive production cross sections for pions hyperons, kaons for proton-proton 
interaction as functions of cm. energy. Solid lines are the results obtained from our model. Data from Ref. p5| . 

In addition, the present model also describes well the energy dependence of the particle production cross sections, 
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as shown in Fig. Here we show the incident energy dependence of the inclusive data for tt"*", tt", tt", K'^^K^.A, 
I]~ and S*^ productions from proton-proton interactions in comparison with the experimental data |25|] . 

The comparisons shown until now in Figs. 2 - 13 show that the combination of particle production mechanisms by 
the resonance decay and the string decay enables us to explore a wide incident energy region from a few hundred MeV 
to a few ten GeV, with reasonably well fitted inclusive as well as exclusive cross sections, which are essential inputs 
in cascade models. 



III. RESULTS 

In the following, we systematically apply our hadronic cascade model (JAMl.O |^^) to proton, silicon and gold 
induced reactions at AGS energies and investigate the effect of cascading by comparing the results obtained by the 
cascade model with the Glauber type calculations. 



A. COMPARISON TO E802 DATA 



In this section, we first focus our attention on the proton transverse distributions to check the detailed examination 
of the collision term and its space-time picture (formation time) used in our model and also to see the transition of 
the reaction dynamics from proton induced collisions to heavy ion collisions. 



proton distributions 
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FIG. 15 

sections 

protons from p+Bc 
p-|-Cu, p-|-Au and 
Si-hAl (6 < 1.79fm), 
(6 < 2.2fm), and 
{b < 2.9fm 
14.6GeV/c. 



Invariant cross 
of 



P+Al, 
central 
Si+Cu 
Si-fAu 
collisions at 
The calcu- 
Cascade 



lated results from 
model (histogram) are com- 
pared with the E802 data 
from Ref. For pro- 

ton induced collisions (up- 
per panel), rapidity interval 
is y = 0.7 (bottom spec- 
trum) to y = 2.3 (top spec- 
trum) with 5y — 0.2. For 
Si-|-Al and Si-|-Cu collisions, 
rapidity interval is y = 0.4 
to y = 2.0 with 5y = 0.2. 
For Si+Au collisions, rapid- 
ity interval is y = 0.4 to 
y = 2.0 with 5y = 0.2. The 
spectra are increased by a 
factor of 10 from bottom to 
upper. 



We show in Fig. ^ proton invariant transverse mass distributions calculated by our cascade model for the pro- 
ton induced reactions, p-t-Be, p-l-Al, p-|-Cu, p-|-Au, and silicon induced reactions, central Si-|-Al, Si-|-Cu, Si-|-Au at 
14.6GeV/c in comparison to the data from the E802 collaboration |39 40]. In each figure, spectra are plotted in a 



rapidity interval of 0.2 and are displayed by multiplying each by a power of 10 from bottom to upper. Si-|-Al, Si-|-Cu 
and Si-|-Au data correspond to the central collision with 7 % centrality. For the calculations of Si-|-A (A=Al,Cu,Au) 
systems, impact parameter are distributed b < 1.797fm for Si-|-Al, b < 2.2fm for Si-|-Cu and b < 2.9fm for Si-t-Au. 
Our calculations show good agreement with the data in proton induced reactions. In silicon induced reactions, our 
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calculations well account for the experimental data in general. However, we can see some overestimate at low trans- 
verse momenta, in particular for Si+Au system. As a result, our cascade model gives larger proton stopping than the 
data. 




Now we compare the cascade model results with the Glauber type calculation in order to see the effect of pion 
rescattering, nucleon cascading, and the consequent deviation from the linear extrapolation of sum of binary nucleon- 
nucleon collisions to proton-nucleus and nucleus- nucleus collisions. The Glauber type models such as the FRITIOF 
model |l] have been widely used at higher energies, i.e. more than 200 AGeV. We use the same method as the 
FRITIOF model with some modifications. The wounded nucleons become resonances or strings in each nucleon- 
nucleon collisions, and strings can interact again before they fragment. Resonances can be converted to nucleons, and 
strings are allowed to be de-excited to the minimum string mass. After all binary collisions are completed, strings and 
resonances are forced to decay. Rescattering of produced particles are not considered. In a present treatment, we have 
used our parameterization in calculating the probability to excite nucleons to resonances or strings. Figure ^ shows 
the results obtained by this Glauber type calculation. In beam rapidity region, for all systems, good agreements can 
be seen because the effect of rescattering would be small as expected. For proton-induced reactions, the Glauber type 
calculation gives steeper shape in comparison to the data at mid-rapidity and target regions. In heavy ion reactions, 
this deviation is significant at around the mid-rapidity. Rescattering, therefore, is necessary to account for transverse 
momentum distributions of protons for reactions involving heavy nuclei. 
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Positive pion distributions 
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FIG. 17. Invariant cross 
sections of positive pions 
from p+Be (0.9 <y< 2.7), 
p+Al (0.7 < y < 2.7), p+Cu 
(0.9 < y < 2.7), p+Au 
(0.7 < y < 2.7), Si+Al 
(0.6 < y < 2.8), Si+Cu 
(0.6 < y < 2.8) and Si+Au 
(0.7 < y < 2.7) collisions 
at 14.6GeV/c. The calcu- 
lated results from cascade 
model (histogram) are com- 
pared with the E802 data 
from Ref. The bin 

of rapidity is 0.2 for all col- 
lision systems. The spectra 
are increased by a factor of 
10 from bottom to upper. 
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Negative pion distributions 
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FIG. 18. Invariant cross 
sections of negative pions 
from p+Be (1.1 <y< 2.7), 
p+Al (0.9 < y < 2.7), p+Cu 
(1-1 < y < 2.7), p-fAu 
(0.7 < y < 2.7), Si-fAl 
(0.6 < y < 2.8), Si+Cu 
(0.6 < y < 2.8) and Si+Au 
(0.7 < y < 2.7) collisions 
at 14.6GeV/c. The calcu- 
lated results from Cascade 
model (histogram) are com- 
pared with the E802 data 
from Ref. The bin 

of rapidity is 0.2 for all col- 
lision systems. The spectra 
are increased by a factor of 
10 from bottom to upper. 
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We now turn to the mass dependence of the pion transverse distributions. In Fig. and Fig. ^ we show the 
calculated and n~ spectra by histograms together with the E802 data |39y4C| ] . Agreement between the cascade 
model and the data is very good for all the combination of projectile and target, and for both of the slope parameter 
and the absolute value of the cross section. 

On the other hand, as shown in Fig. |l^ and Fig. Glauber type calculations well reproduce the data for small 
mass systems p+Be, p+Al, p+Cu and Si+Al, and give similar slopes to the experimental data, while the multiplicity 
of pions in heavy systems are larger than those of the data. This is due to the effect of rescattering in which pious 
are absorbed during the evolution for the large mass systems. 



Positive pion distributions (Glauber type calculation) 
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FIG. 19. Same as 

Fig. |l^, but the histograms 
show the results obtained 
from the Glauber type cal- 
culations with- 
out any rescattering among 
produced particles. 
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Negative pion distributions (Glauber type calculation) 
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FIG. 20. Same as 

Fig. |l^, but the histograms 
show the results obtained 
from the Glauber type cal- 
culations with- 
out any rescattering among 
produced particles. 



Positive kaon distributions 
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FIG. 21. Compari- 
son of the invariant trans- 
verse momentum spectra of 
K'^ between the cascade 
model and the experimen- 
tal data of 11,1^ in p-hBe 
(0.9 <y<2.l with bin size 
0.2), p-hAl (0.9 < y < 2.1 
with bin size 0.2), p-fCu 
(0.9 < y < 2.1 with bin size 
0.2), p-fAu (0.7 < y < 2.1 
with bin size 0.2), Si-|-Al 
(0.6 < y < 2.2 with bin size 
0.2), Si-HCu (0.6 < y < 1.8 
with bin size 0.4) and Si-|-Au 
(0.7 < y < 2.1 with bin size 
0.2) collisions at 14.6GeV/c. 
The spectra are increased by 
a factor of 10 from bottom to 
upper. 
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Negative kaon distributions 




FIG. 22. Comparison 
of the invariant transverse 
momentum spectra of K~ 
in p+Be (1.1 < j/ < 1.9 
with bin size 0.2), p+Al 
(0.9 < y < 2.1 with bin size 
0.2), p+Cu (0.9 < y < 2.1 
with bin size 0.2), p+Au 
(0.9 < 2/ < 1.9 with bin size 
0.2), Si+Al (1.0 < 2/ < 1.8 
with bin size 0.4), Si+Cu 
(1-0 < y < 1.8 with 
bin size 0.4) and Si+Au 
(0.9 < y < 2.1 with bin size 
0.4) coUisions at 14.6GeV/c. 
The spectra are increased by 
a factor of 10 from bottom to 
upper. 
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Positive kaon distributions (Glauber type calculation) 




FIG. 23. Glauber 
type calculations of invari- 
ant cross sections of for 
p+Be, p+Al, p+Cu, p+Au, 
Si+Al, Si+Cu and Si+Au 
reactions at 14.6GeV/c. in 
comparison to the E802 data 
from 11,0. The mean- 
ing of figure is the same as 
Fig. 
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Negative kaon distributions (Glauber type calculation) 




FIG. 24. Glauber 
type calculations of invari- 
ant cross sections of K~ for 
p-l-Be, p-l-Al, p-l-Cu, p+Au, 
Si-I-Al, Si-I-Cu and Si+Au 
reactions at 14.6GeV/c. in 
comparison to the E802 data 
from The mean- 

ing of figure is the same as 
Fig. H 
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Let us study the kaon and anti-kaon transverse mass spectra of E802. The calculated transverse mass distributions 
of and K~ in cascade model are compared with the E802 data [^,^ in Figs. ^ and Figures ^ and |2^ 
are the results with the Glauber type calculations for and invariant transverse momentum distributions. We 
find more significant differences between the cascade and Glauber results of kaon pr oduc tions than those in proton 
and pion spectra. This fact shows the importance of the rescattering: As discussed in [ID, some of the exoergic AIB 
reactions, which involve resonances, have very large strangeness production cross sections, and they contribute to 
K and K productions, especially in heavy ion reactions [T^j3l[ |. Enhancements due to these meson rescattering are 
clearly seen in if + and K~ spectra, except for the reactions of p-|-Be and p-f Al, because there is no meson-baryon 
collision in Glauber type calculation. 
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B. Au+Au COLLISIONS 



Au+'^'Au^p+x atll.6AGeV/c 




FIG. 25. The result 
of Cascade model calcula- 
tion of transverse mass dis- 
tributions of protons for 
central Au-|-Au collision at 
11.6GeV/c in different ra- 
pidity intervals. The spec- 
tra are scaled down by a fac- 
tor of 10 successively from 
upper corresponds to the 
cm. rapidity y — 0.05 to 
lower spectrum of y = 1.05 
with the bin width of 0.1. 
Impact parameter is dis- 
tributed from to 3.338fm. 
Data are from Ref. pl[ . 



We continue our comparison to E866 experimental data with the truly heavy ion collision Au-|-Au at ll.GAGeV. 
Our cascade model calculation with impact parameter b < 3.338fm is compared to E866 data |^ in Fig. ^ from cm. 
rapidity of y = 0.05 to y = 1.05 with the rapidity bin of 0.1 scaled down by a factor of 10 from the top to the bottom 
spectrum. The cascade model results show good agreement with the data at the transverse mass above 0.2GeV/c^. 
The cascade model, however, overpredicts again the proton spectrum in the low transverse momentum region. 
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FIG. 26. Comparison of rapidity distributions of protons 
(circles), positive pious (squares) and negative pious (tri- 
angles) between experimental data and cascade model 
calculations (histograms) for central Au-|-Au collision at 
11.6GeV/c. The data for positive pions is scaled down by a 
factor of 0.5 and for negative pions 0.25. In cascade calcula- 
tions, impact parameter is distributed from to 3.338fm. 
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FIG. 27. Transverse mass distributions of protons and 
pions for central Au-|-Au collision at 11.6GeV/c. The trian- 
gles represent the data of protons scaled by a factor of 0.1 
squares correspond to ■k~ data and circles data scaled 
by a factor of 10 from from Ref. ^ij. Histograms repre- 
sent the results from Cascade model with impact parameter 
6 < 3.338fm. 
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In Figs. and [2^, we compare the cascade model results with the experimental data by E802 collaboration |Q in 
central Au+Au collisions. Pion multiplicities are in good agreement with data as well as the slopes of both tt"*" and 7r~ . 
However, the present cascade model does not describe the suppression of protons having low transverse momenta, 
and consequently it gives stronger stopping of proton than the data. This proton rapidity spectrum for central 
Au+Au collisions shows a similar amount of stopping as those with other cascade models like RQMD ARC |T^ ] 
and ART |p^ . In addition, Glauber type calculation gives the same results for the proton rapidity distribution. 
Therefore, this defect is not a consequence of the cascade model. Since the deviation of the transverse mass spectrum 
of heavy hadrons from a single exponential behavior is generally considered as a result of the radial flow, it may be 
influenced by the nuclear mean field. In fact, it has been found in the works of RQMD [Q and ART nuclear 
mean field significantly reduces maximum baryon densities of the hadronic matter, and consequently the midrapidity 
protons becomes small. 

In this work, we have assumed that the elastic and inelastic cross sections involving baryon resonances, except for 
the de-excitation processes to NN, are the same as that in NN channel at the same cm. momentum. However, since 



the de-excitation cross sections are enhanced due to the generalized detailed balance as explained in Section |IIC| , if 
other cross sections is smaller than those of NN, stopping power may be reduced. 



C. Mass dependence of the collision dynamics 

In this section, we study the mass dependence of the collision dynamics within the hadronic cascade model. First 
we present the cascade model results of BB, MB and MM collision number as a function of system mass {A + B) in 
comparison to that of Glauber type calculations. 




mass number mass number (A+B) 



FIG. 28. Mass dependence of number of BB (open squares), MB (open circles) and MM (open crosses) collisions obtained 
from cascade model calculations and Glauber type calculations (open diamonds) for p+A collisions in the left panel and for 
Si-|-A and Au-(-Au collisions (right panel). 

Figure ^ displays the total collision number of BB (open squares) and MB (open circles) and MM (open crosses) 
obtained from cascade calculation together with Glauber type calculation (open diamonds) for the p-l-A (left panel) 
and Si-|-A (right panel) collisions. When system becomes bigger, BB collisions are much more frequent in cascade 
model than in the Glauber predictions even in the proton induced collisions. This indicates that there are successive 
nucleon cascading in the nuclear medium in the cascade model picture. It is interesting to see that the number of 
BB and MB collisions are almost the same in heavy ion collisions. This seems to be the origin of the pion number 
suppression, the increase in proton transverse momentum slope, and the increase in kaon yield. 
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FIG. 29. Mass dependence of number of total mesons from cascade model calculations (full circles) and Glauber type 
calculations (full triangles) for Si+A (A=A1, Cu, Au), Au+Au collisions (right panel) and for p+A (A=Be, Al, Cu, An) 
collisions (left panel). Meson multiplicity is reduced by including rescattering. 

Indeed we can see the reduction of produced pion multiplicity in cascade model compared to the Glauber type 
calculations as shown in Fig. 29 where the number of total produced mesons are plotted as a function of system mass 
number. Pions are absorbed mainly in the two-step processes in the cascade model. For example, the most important 
pion absorption path at AGS energies is 

TriV -> A , iVA NN . (39) 
Therefore, large number of BB and MB collisions are necessary to describe appropriate pion absorptions. 




FIG. 30. Collision spectrum of BB (full lines), MB (dashed lines) and MM (dotted lines) collisions from the cascade model 
calculations and Glauber type calculations (crosses) for the p+Be, p+Al, p+Cu, p+Au, Si+Al, Si+Cu, Si+Au and Au+Au 
collisions at AGS energies. 

In order to get more detailed information on the mass dependence of the collision dynamics at AGS energies, we 
display in Fig. |3^ the colliding energy spectrum of BB, MB and MM given by the cascade model together with 
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those from the Glauber typer calculations. The BB collision distributions as a function of invariant mass are very 
different between the cascade model and the Glauber type model. The Glauber type calculations predict the collisions 
which are spread around the initial NN cm. energy, while the BB collisions occur at all available collision energies 
in the cascade model. In both p+A and A+B reactions, the MB collisions are pronounced in the resonance region 
(•^i < 2GeV). It is interesting to note that in p+Au system, the number of low energy BB collision is much larger 
than that of p+Be, p+Al and p+Cu systems. In A+B systems (bottom of Fig. |3^), collision number grows very 
quickly, however, the shape of the collision spectrum is similar in all the systems. 



IV. SUMMARY 



We have systematically studied the system mass dependence of the particle distributions at AGS energies with 
a newly developed cascade model (JAMl.O). The cascade model is shown to provide a good description of the 
observed data for various combination of projectile and target without any change of model parameters. The effect 
of rescattering of produced particles and nucleon cascading are found to be important to explain both the pion yield 
and the transverse slope, which are demonstrated by comparing the cascade model results with the Glauber type 
calculations. Those effects increase the transverse momentum slopes of protons and pions, and reduce the pion yield. 
The importance of the rescattering among particles is more visible in kaon spectra. 

One of the problems in the hadronic cascade model JAM is that it gives much larger stopping of the protons 
in nucleus-nucleus collisions. This large stopping is not the consequence of the rescattering because Glauber type 
calculation also gives the same amount of baryon stopping. This problem of strong baryon stopping in cascade models 
has been reported that proton spectra can be fitted by the inclusion of nuclear mean field in RQMD and ART . 

Another possible solution may be to suppress the cross sections such as cr(iVj*A^| N^N^), which is assumed to be 
the same as that in NN ingoing channel at the same cm. momentum in a present model. These interactions among 
resonances become important at AGS energies where we have sufficiently dense matter in heavy ion collisions, and the 
baryon stopping is sensitive to the cross sections in the resonance ingoing reactions. In fact, we have checked that if 
resonance-resonance (BB) cross sections are reduced from nucleon-nucleon cross section, we get less proton stopping 
than the present results. Detailed study in this line will be interesting. 
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